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We study a Hamiltonian lattice version of the two-dimensional Wess-Zumino model. Preliminary results ob- 
tained by Quantum Monte Carlo with a many-parameter guiding wave function are presented. We analyze the 
pattern of supersymmetry breaking by measuring the ground state energy and a set of supersymmetric Ward 
identities. The algorithm is quite effective and allows very precise measurements. 



1. INTRODUCTION 

Numerical simulations of lattice field theories 
are usually performed in the Lagrangian formu- 
lation. The alternative approach of the Hamil- 
tonian formulation ^ has received very little at- 
tention so far. Nonetheless, we think there are 
very good reasons to develop numerical simula- 
tion techniques for the Hamiltonian approach: 
powerful many-body techniques are available 
which allow the direct computation of the vac- 
uum wave function properties; fermions are im- 
plemented directly and need not be integrated 
out; property like the mass spectrum are more 
immediate. A check of universality between the 
Lagrangian and the Hamiltonian formalism would 
also be very welcome. 

2. THE MODEL 

The continuum two-dimensional Wess-Zumino 
model is defined by the supersymmetry charges 
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where 4'{x) is a real scalar field, p{x) is its conju- 
gate momentum, and ipi^) is a Majorana fermion. 



The charges obey the = 1 supersymmetry al- 
gebra 

Ql = Ql = P" = H, {Qi, Q2} = = 2P, 

which we have written in the Majorana basis 
70 = C = —(72, 71 = ica. Since P is not con- 
served on the lattice, a lattice formulation of a 
supersymmetric model must break the algebra 
explicitly. A very important advantage of the 
Hamiltonian formulation is that, since H is con- 
served exactly, we can maintain a 1-dimensional 
supersymmetry subalgebra, e.g., Q\ = H, and 
expect to recover the rest in the continuum limit, 
together with Lorentz invariance |^; this subal- 
gebra is enough to guarantee some of the most 
important property of supersymmetry, including 
pairing of fermionic and bosonic states of nonzero 
energy; spontaneous breaking of supersymmetry 
is equivalent to nonzero ground-state energy Eq. 
We adopt the lattice formulation Q 
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with canonical (anti)commutation rules, and H ~ 
. We can replace the Majorana fields V'^'^ by 
one Dirac field x by the clever transformation Q 
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obtaining the Hamiltonian 

L 



Pn 
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v{ci>n: 



(xJ.Xn+1 + h.c.) + (^„) {2xlxn - 1) 



which is expHcitly free of sign problems, except 
potentially at the boundary; we adopt free bound- 
ary conditions to avoid the sign problem alto- 
gether n. 

In strong coupling, the model reduces to a su- 
persymmetric quantum mechanics for each site; 
supersymmetry is broken if and only if the degree 
of the prepotential V is even . In the continuum 
(and on the lattice in weak coupling), supersym- 
metry is broken at tree level if and only if V has 
no zeroes. The predictions of strong coupling and 
weak coupling are quite different, and it is inter- 
esting to study both numerically and analytically 
the crossover from strong to weak coupling. 

As a benchmark of the case where supersymme- 
try is unbroken, we consider the cubic prepoten- 
tial V = (f>^. In the more interesting case of bro- 
ken supersymmetry we study the quadratic pre- 
potential V = (p+Xo- Here, spontaneous symme- 
try breaking is expected to occur at Aq > 0, while 
for Ao < the scenario is less clear; in fact, weak 
coupling perturbation theory predicts a symmet- 
ric ground state while strong coupling predicts 
broken supersymmetry with Eq ~ exp(cAo) for 
Ao — oo. 

3. MONTE CARLO SIMULATIONS 

We perform our simulations by the Green Func- 
tion Monte Carlo algorithm with the guiding 
wave function associated to the trial ground state 

|*o)trial =exp J^n^-ri |^'o)froo, 

Fn=Pl{<Pn) + i-irixiXn - l/2)P2(0„), 

where |^'o)free is the free {V = 0) ground state 
and P\^2 are polynomials optimized according to 
the algorithm described in Ref . ; see Ref. [| for 
more details. We wish to remark that, in order 
to keep the variance of observables finite as the 
simulation proceeds, it is necessary to simulate a 
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Figure 1. Eq/L vs. l/K for V 
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population of K walkers (field configurations at 
fixed time), and extrapolate the results to K ^ 
oo. 

The trial wave function is an important fea- 
ture of our simulations, and we are able to de- 
termine accurately the free parameters. The trial 
wave function allows a substantially improvement 
of the quality of our numerical results; moreover, 
since it is an approximation to the exact ground 
state, it contains important physical information 
about the model. 

This algorithm can be parallelized effectively 
on a network of PCs connected through Ethernet; 
our MPI code reaches 90% efficiency. 

We measure the ground-state energy i?o, which 
is the simplest and most precise observable to 
measure in a Quantum Monte Carlo, and a num- 
ber of observables of the form {Q,Xq}, Xg = 
J2n '/'nV'nj which are zero for unbroken supersym- 
metry (for any fermionic operator X, {{Q,X}) = 
is a supersymmetry Ward identity). 

Fig. d shows the ground-state energy per site 
as a function of l/K for the cubic prepotential 
V — (j)'^ and L — 10. The evidence for unbroken 
supersymmetry is quite convincing. The bosonic 
and fermionic contribution to Eq/L are ~ ±0.7: 
we are observing a cancellation of the order of 
lO^'*. Similar conclusions can be drawn by look- 
ing at {Q,Xq}, with smaller numerical accuracy. 



Fig. H shows the ground-state energy per site 
as a function of 1 / K for the quadratic prepoten- 
tial = 02 + Ao, for Aq = -1, -1.25, -1.5, 
and —2 and L = 10, 22. The data are insensi- 
tive to the value of L, and they are consistent 
with a hnear extrapolation to \/K 0. Super- 
symmetry breaking is clear for Ao > —1.25; for 
smaller values of Aq, we cannot discriminate be- 
tween £^0=0 and an exponentially decreasing 
Eq. As in the cubic case, the study of {Q,Xq} 
gives similar conclusions, with smaller numerical 
accuracy. 

4. CONCLUSIONS 

The Hamiltonian approach is a worthwhile al- 
ternative to the usual Lagrangian approach to 
lattice field theory. Accurate Monte Carlo sim- 
ulations can be performed and compared to ana- 
lytical predictions, e.g., strong coupling computa- 
tions. Fermions can be implemented directly and 
need not be integrated out, thus avoiding a non- 
local determinant; for certain models in 1 + 1 di- 
mensions, like the one studied here, they present 
no additional costs compared to bosons. It should 
be observed however, that in higher dimension 
fermions lead to sign problems which need to be 
solved 1^. Another important advantage of the 
Hamiltonian formalism is the possibility of pre- 
serving exactly a 1-dimensional supersymmetry 
algebra. 

In the light of the encouraging preliminary 
results presented here, we plan to extend our 
study in several directions. We are computing 
the strong-coupling expansion, and we will study 
the weak-coupling expansion and the perturba- 
tive continuum limit. We will investigate numeri- 
cally the transition from strong to weak coupling, 
looking at the pattern of supersymmetry break- 
ing. 
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